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Abstract 

A non-supersymmetric model of two-loop radiative seesaw is proposed. The model contains, 
in addition to the standard model (SM) Higgs boson, an inert SU{2)i doublet scalar rj and two 
inert singlet scalars (p and x- The lepton number is softly broken by a dimension two operator, 
and the tree- level Dirac mass is forbidden by Z2 x Z2 (or D2n) , which predicts the existence of 
two or three dark matter particles. The scalar sector is minimal; none of the scalar fields can be 
suppressed for the radiative seesaw mechanism to work. There are by-products: The SM Higgs 
boson decay into two 7's is slightly enhanced by r/+ (the charged component of r/) circulating in 
one-loop diagrams for h — t- 77. The 135 GeV 7-ray line observed at the Fermi LAT can be also 
explained by the annihilation of x dark matter. We employ a mechanism of temperature-dependent 
annihilation cross section to suppress the continuum 7-rays and the production of anti-protons. 
The explanation can survive even down to the XENONIT sensitivity limit. 
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I. INTRODUCTION 



Why the neutrino masses are small is a long-standing mystery. The seesaw mechanism 
is an approach to provide an answer to it. The traditional seesaw mechanism indicates 
the existence of a superhigh scale physics beyond the standard model fSM). Another way to 
confront this problem is to generate the neutrino masses radiatively ji-Q]. Many models have 
been proposed on the basis of the radiative generation of the neutrino masses, and the idea of 
the radiative seesaw mechanism 0, [sl is along the line of this idea: Right-handed neutrinos 
are introduced, but the Dirac masses are forbidden by a discrete symmetry. This discrete 



symmetry can be an origin of stable dark matter (DM) particles [7Kll|] in the universe. 

To produce the neutrino masses, the lepton number L has to be violated. In most of 
the models the lepton number L is violated softly by dimension three operators, Majorana 
masses or scalar tri-linear couplings, and the number of the loops is ranging from one to 



three 



see 



12 



for different models). For the radiative seesaw mechanism, the number of 



loops i means a scaling down of (1 -|- 3i) orders of magnitude for the right-handed neutrino 
mass (see also the discussions of [l3|); {n/lGn'^Y ~ (k/0.1) x 6.3 x 10"'' , (k/0.1)^ x 4.0 x 
10~^ , (k/0.1)^ X 2.5 X 10~^° for i = 1 , 2 and 3, respectively, where k is a generic coupling. 
Since the Majorana mass of the right-handed neutrino for the tree-level seesaw is 0(10^") 
GeV, we may obtain a Majorana mass of 0(1) TeV naturally in two-loop radiative seesaw 
models 

In this paper we propose a radiative seesaw model, in which the lepton number is softly 
broken by a dimension two operator, and the neutrino masses are generated at the two-loop 
level. The discrete symmetry is Z2 x Z2 (or D2N with = 2,3, . . . '). Therefore, two or 
three DM particles can exist in this model ^. Obviously, radiative generation of the neutrino 
masses means an extension of the SM Higgs sector, which may have impacts on the existing 
experiments. In our model, we have a set of an inert doublet scalar rj, and two singlet scalars. 
This set is minimal in the sense that the radiative neutrino mass generation does not work 
if one of them is suppressed. So, none of the extra scalar fields is ad hoc introduced. The 
existence of additional scalar doublets can change the decay rates of the SM Higgs boson 
h. The results of LHC indicate a slight excess of /i — )■ 77 20l-l22| which in fact could be 
explained by an additional inert doublet circulating in one-loop |23l-l26|. 

There is yet another excess of 7 at the Fermi Large Area Telescope (LAT) 
There are analyses 



27-30 



31 



36 



that indicate a monochromatic 7-ray line of 135 GeV in the 
Fermi data. It has been reported 37| that a two-component DM system consisting of an 
inert doublet scalar and a scalar can explain the monochromatic 7-ray line at Fermi LAT. 



Different models with a two-component DM have been also considered in 38|, [SOj to explain 



^ D2N (the dihedral group of order 2N) is larger than Z2 x Zt^. However, we use only the one-dimensional 
representations of D2N so that D2N acts as Z2 x 



Multicomponent DM system has been considered recently in 16l4l9l|. and see the references therein. 
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TABLE I: The matter content of the model and the corresponding quantum numbers. Z2 x Z'^ is the 
unbroken discrete symmetry, while the lepton number L is softly broken by the ^ mass. Din (N = 2,3,...) 
is the dihedral group of order 2N. 
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the monochromatic Fermi LAT 7-ray hne. Therefore it is natural to wonder whether our 
two-loop radiative seesaw model can explain the 7 excess in the Higgs boson decay as well 
as at the Fermi LAT ^. We find that this is in fact possible if we accept that certain scalar 
couplings are large at the border of perturbation theory, where to suppress sufficiently the 
continuum 7's and the production of anti-protons, we employ a mechanism of temperature- 
dependent annihilation cross section 

II. THE MODEL 

The matter content of the model is shown in Table I. The new fields are (in addition 
to the right-handed neutrino N"^): SU{2)l doublet scalar r] (L= —1), and singlet scalars x 
(L= 0) and (L= 1), where L is the lepton number. Note that L of is zero and that 
four different representations of Z2 x Z2 are exactly the singlets of the dihedral group of 
order 2N, D2N {N = 2, 3, . . .). The Z2 x Z2 x L -invariant (or D2N x L -invariant) Yukawa 
couplings of the lepton sector can be described by 

Cy = Y,'^,H^L,l'nj + Y,lL,evN^^, + h.c. (1) 

with the Majorana mass term of the right-handed neutrinos iV|.^ {k = 1, 2, 3) 

CM.i = -l[M,N^^,N^^, + h.c.] . (2) 
The most general form of the Z2 x L-invariant scalar potential is given by 
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Recently, it has been argued that the monochromatic 7-ray line can be explained basically by the same 



one- loop contribution as for — > 77 [37|, |40|443 [ 
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+7ix' + l2{H^H)x^ + 73(r?^r7)x' + 74I0I' + l,{H^Hm^ 
+l,{v^v)\<P? + l-rx'\<P? + ^[{H^V)x^ + h.c.] . (3) 

Note that the "A5 term", {l/2)X^{H''ri)^, is forbidden by L. The Z2 x Z2- invariant mass 
term is 

Vm = mlH^H + mlr]^r] + ^mlx' + ml\(P\^ + ^ml[^^ + {ff] , (4) 

where the last term in (jlj) breaks L softly. This is the only Z2 x -invariant mass term 
which can break L softly. In the absence of this term, there will be no neutrino mass. The 
charged, CP even and odd scalars are defined as 

The tree-level masses of the scalars are given by 

= 2Aiw^ , mj± =ml + ^Xsvl , mjo^ = mjo = + ^(As + X4)vl , (6) 
2 2,2, 2 2 2 2, 2 2 2, 2 



A. Stability of the vacuum 

The potential Vx is unbounded below if 



Ai , A2 , 71 , 74 > , (7) 

A3 > -^\/a^,A3 + A4 > -^V^A^, (8) 

2 / 2 / 2 / 

72>-gV^i^i' 75>-2V^i^4, 73>-2V^27i, (9) 

76 > -gV^Ml , 77 > --VTiTI , (10) 



2 



Ai + A2 + 7i + 74 - 2 ( V A1A2 + V Ai7i + \l A174 



+ V A271 + V A274 + A/7i74j > \tA (11) 

are satisfied. The minimum of V\ is zero if the inequalities above are satisfied. The dis- 
crete symmetry Z2 x Z'2 (or D2n) is unbroken if, in addition to fl7|)- f|TT]) . the inequahties 
m\ , m\ , m| and m\ — \m^\'^ > are satisfied. 

The inequalities (l7|)-f lTT]) are sufficient conditions, but not necessary ones. If we assume 
that ([7j) and 72,73,75,76,77 > are satisfied, for instance, the inequalities (IHj) are relaxed 
to 



A3 , A3 + A4 > -2a/AiA2 , (12) 
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and (fTT]) is relaxed to 

Ai + A2 + 7i + E7*-2\/a^> |/«| . (13) 

i=2 

Since nih = 125 GeV and Vh = 246 GeV, the Higgs coupling Ai is fixed at 0.129. Then f fT2|) 
implies that 



As , A3 + A4 > -2.5v/A2/47r. (14) 



B. Neutrino mass 



The neutrino masses can be generated at the two-loop level as shown in Fig. [T] The 
mechanism of the radiative generation is the following. Because of the soft breaking of the 
dimension two operator 0^, the propagator between and can exist. This can generate a 
^0^0 ]2iass term. In the one-loop radiative seesaw model of this mass is generated at the 
tree-level through the "A5" coupling. So the effective \^ is 

(15) 



A 



eff ^ ____ 

^ 647r2 



2 2 2 2 

m, m, m, m, 

VI 1„ 91 VR. 1 VR 

2 9 9 2 9 9 



Therefore, Ag*^ can not be large within the framework of perturbation theory (A|^ ^ 
0(10^'^)). The two-loop neutrino mass matrix is calculated to be 

{My)ij = [jq^) ~^J2'>^ikY/kMk dx{ Bo{-x,m^,m^^) - Bo{-x,m^,m^^) } 

X 

TTT for m„ = m„o ~ m„o (16) 



(x + m2)2(x + M|) 

yi^yi^ / /TTi \ 2 \ 

E iM 1^ ) + « ' (17) 



^ 167r2Mfc \ V M, 
where the function Bq is the Passarino-Veltman function [45 1 

-^5o(p^ml,m2) = / 7^-^77^ 2 I • N^/;^ I — ^2 2 , • a ' (1^) 

and to obtain (fT7|) we have used the one-loop mass formula with the effective couphng ( 1TB]) 
^. Therefore, the scale of the light neutrino mass will be 

2 1 2 ^ s 9 2 

K 1 mj) f K Y -7 



10-7 X ^ fl9l 
647r2l67r2 M VO.l/ M ' ^ ^ 

where rajj/M is the scale in the case of the tree-level Type-I seesaw. This means that we 
can scale down the mass of the right-handed neutrino by seven orders of magnitude. So, the 
right-handed neutrino masses of TeV or less are naturally expected in this model. 



There is 0(1) correction to the approximate formula (|17p which we have checked numerically. 
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FIG. 1: Two-loop radiative neutrino mass. 



C. Constraints 



1: — )■ e 7 

The constraint coming from /x — )• 67 is given by 

3a 



46 



B{fi 67) 



2 ^2 



k 



'Ml 



6(1 -x) 



'1 - 6a; + 3x^ + 2x^ - 6x^ In x) 



< 2.4 X 10 



12 



(20) 



where the upper bound is taken from 47|]. A similar, but shghtly weaker bound for 
r — 7- /i(e)7 given in [47| has to be satisfied, too. Since -^2(2;) ~ l/Sx for a; ^ 1, while 
1/12 < F2{x) < 1/6 for < X < 1, the constraint can be readily satisfied if Mk ^ or 
Mfc > m^±. If we assume that Mi = M2 = M3 = 1 TeV > m^± in ([201), the constraint (I2U]1 
becomes 5(/i ^ 67) ~ 10"^ x Efc ^^^^^fc ^ 2.4 x IQ-^^^ Therefore, Y^IY/^^ < 0(10"^) can 
satisfy the constraint. This size is also the right size to obtain light neutrino masses of a 
correct order. 



2: g,~2 

The extra contribution to the anomalous magnetic moment of the muon, = {g^ — 2)/2, 
is given by 46 1 



5a, 



If we assume that | ^atXHk^'^ 



m, 



167r^m^_i_ 

Mi 



'Ml 



Mi 



(21) 



I, then we obtain 



2.2 X 10~^5(/i ^ 67) ^ 3.4 X 10^ 



^11 



(22) 
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FIG. 2: One-loop diagrams for h — > 77. 



if f l2U]) is satisfied, where the upper bound is taken from 48|]. So, the constraint from 
has no significant influence. 



3: Electroweak precision 

The electroweak precision measurement requires 



48 



49 



AT ~ 0.54 



0.02 



-0.11 
-0.12 



(23) 



for rrih = 115 - 127 GeV. Therefore, |mj,± — "^,,0 1 , |m^± — m^o| ^ 100 GeV is sufficient to 
meet the requirement. 



III. DM AND RESTRICTING THE PARAMETER SPACE 



A. h 



77 



Because of the coupling Cjj^jj^^ = —\^{WH){rj^vi) = —X3Vhhri~^r]~ + ■ ■ ■, there are 
additional diagrams that contribute to the one-loop decay h — )■ 77,^ which are shown in 
Fig. [2J Therefore, the decay width for two 7's can be increased 



23 
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r(77) 

FSM(77) 

n 

Fi/2{t) 

Fiir) 
Foir) 



3(3/2)^Fi/2(r,) + F^jrw) + 2A3(m^/(7^mg±)Fo(r^±) 
3(3/2)2Fi/2(n) + Fi(rH.) 
Ami/ ml, Tw = 4:mly/ml , r^± = Am\±lm\ , 
2 + 3r + 3r(2 - r) arcsin2( 1/^/7) , 



-2ril 



r arcsm 



:i/v^)] 



"[1 — r arcsin^(l/^/r)] 



(24) 
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[ GeV ] 



FIG. 3: The area with r/F^'^ — 1.6 ± 0.4 in the m^i-As plane. The horizontal red line is the stability 
bound in given in (|14p . 

where arcsin^(l/Y^) should be replaced by (— l/4)[ln — ztt]^ for r < 1. We obtain 



r(77) 



1.20 1.31 1.54 

1.17 1.26 1.45 

1.16 1.24 1.42 

1.05 1.12 1.22 



for A3 = — { 1 1.5 2.5 and m^-t = < 



135 
145 
150 
200 



GeV. (25) 



So, if the charged inert scalar 77"*" is relatively light and A3 is negative and large, the observed 
excess 1.6 ± 0.4 in the CMS experiment 
for this mode in the ATLAS experiment 



21 



20 



can be explained. (The best fit signal strength 
2^ is /i = 1.9 ± 0.5.) In Fig. [3] we show the 
area in the m.^±-A3 plane in which F(77)/F^'^(77) = 1.6 ± 0.4 can be obtained ^ As we can 
see from (jS]), a large negative A3 may endanger the vacuum stability, because Ai is fixed 
at 0.129. Therefore, we assume that all the quartic scalar couplings except A3 and A4 are 
positive and use the second set of the inequality conditions ( 1T2|) and (|T3|) . (|T4|) means that 
A3 ^ —2.5 if A2 is at the border of perturbation theory. Thus, the Higgs boson decay mode 
/i — )■ 77 prefers the parameter space: 

1. All the quartic scalar couplings except A3 and A4 are positive and A2 is large, and 

2. the constraints 



- 0.2 < A3/47r < -0.08 and m^± < 200 GeV 



(26) 



The upper bound m^± < 135 GeV given in [26'] is obtained from r{-f-f)/T^^ > 1.3 and Xs/An > 
-1.46/47r ~ -0.116, which is consistent with Note that we use 2.0 > T{"f"f)/r^^{j-/) > 1.2 for 

Fig. [3] along with the stability constraint (fH)) . 
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FIG. 4: Tree-level DM annihilations into W+W^ and ZZ. 



have to be satisfied. 



B. Direct Detection of DM 



As we can see from Table I, eitlier N'^ or 77 can be a DM candidate. Here we assume tfiat 
r]^, tlie CP even component of r], is a DM particle and assume that Mk ^ m^o^ , m^o , m^± 
to satisfy the yU — 67 constraint (1201) . The model can have three stable DM particles in 
principle, but to simplify the situation we assume a two-component DM system. Another 
one is either x 0- As we see from the potential ([3]), there is no significant difference 
between x ^"^^ ^-s DM. So we assume here that x is the second DM particle. 



The spin-independent elastic cross section off the nucleon cr(x(?7^)) is given by 



49 



1 ( 72(AL/2)/mjv \ / "^7v"^x«) 



TT 



(27) 



where A r = A3 + A4, m-Ar is the nucleon mass, and / ~ 0.3 stems from the nucleonic matrix 



element 



50 



The cross sections have to satisfy 



where Qrh^ = 0.116 (5l| and auB(m) ~ 3 x 10"^^ cm^ [52| is the XENONIOO limit for the 




(28) 



DM mass of 135 GeV. So, we find 



I72I ^ 0.035 (0.003) if n^h'^ ~ rirh'^ = 0.116 and = 135 GeV 



IA3 + A4I ^ 0.069 (0.006) if n^jh'^ ~ rirh'^ 



0.116 and m„o 

'IR 



135 GeV 



(29) 
(30) 



where we would obtain the numbers in the parentheses when the XENONIT sensitivity 
a ~ 2 X 10"^^ 



cm 



53( 1 would has been reached. Because of (l26l) we have to assume a 
large negative A3, and, because of m^± — m^o = —\iv\/2, A4 has to be negative, too, to 
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ensure m^± > m^p^. If A4 is negative and in addition A3 is large and negative, the inequality 
( 1T2|) may be violated, unless A4 is small and negative. This means Al — A3, implying that 
( 130|) can not be satisfied. Therefore, the relic density of 77^ has to be small to satisfy the 
constraint This is welcome, because the annihilation cross section of 1]^ is large in 

general due to the gauge interactions shown in Fig. H] Thus, the parameter space is 

further constrained to: 

1. A4 has to be negative and small to ensure m^± > m^o^ and to satisfy the constraint 
IT 



2. since A3 is assumed to be large and negative (see fl26l) ). the constraint fl28l) can be 
satisfied only if Q^i/^t ^ 1, and 

3. to satisfy the XENONIOO (IT) constraint we have to impose I72I ^ 0.035(0.003) for 

= 135 GeV. 



C. Relic densities of DM 

Since our parameter space has been already constrained to a certain amount, we next 
calculate the relic density of DM, Qj- = + Q^i- To simplify the situation we have been 
assuming throughout that rj^ and x a^^e DM particles. In this two-component DM system 
there are three different thermally averaged cross sections 



^0 . 



(31) 



that are relevant for calculating the DM relic density, where SM stands for the SM particles ^. 
No semiannihilation rj^x ^ 0i?(0/)SM is allowed if fn^j^^^^ > m^o +m^. Then the evolution 
equation for Y, the number density over the entropy density, can be written as 

-/iMpL 



16 



dY 

= -0.264 gl/' 
dx 



x^ 



dY^ 



(^(VrVr, XX)v)\Y^^>^Y^o^ 



dx 



-0.264 <^y2 ^ {(a(xx;SM).;)(r,r,-F,F, 



Y Y 



Y Y "^ R '^ R 
-'x-'x / 



(32) 



(33) 



where F is F in equilibrium, x = fi/T, = 1/m^o + l/m^, and T, MpL and = 90 
are the temperature, the Planck mass and the total number of effective degrees of freedom, 
respectively. 



^ These thermally averaged cross sections are tree- level ones and do not include those into two 7's. Anni- 
hilations into two 7's will be separately calculated later on. 



10 



Before we solve the evolution equations numerically, we consider what we would expect. 
As noticed the relic density of 77^ will be very small because of large A3 and gauge interactions 
(i.e. large {(T{r]fjri^; SM)f )), while the annihilation cross section of x into the SM particles 
is suppressed because of f l2^ (i.e. small (cr(xx; SM)t;)). That is, the DM conversion cross 
section (o'(?7^?7^; xx)'^) and the mass difference Am^^ = m^o — will play an important 
role. Note that the smaller ^rrirjx larger is the effect of the DM conversion on Q^. To 

see this more explicitly, we assume that 77^ annihilates very fast so that before and at the 
decoupling of x the 77^ DM is in thermal equilibrium. Then the expression in { } in the rhs 
of (15^ can be written as 



, (34) 



which also appears in the co-annihilation of DM with an unstable particle 0] ^. If — 
m^o < 0, the effective annihilation cross section of x is small at low temperature (large x), 
while it is large at high temperature (small x). Because of the nontrivial interplay between 
jy, and Arrirj^, it may be possible to obtain a correct relic density fi^/i^ = 0.1157 ± 0.0023 



5l| . In Fig. |5]we show the effective annihilation cross section (the expression in [ ] of 
as a function of x = /i/T for m^o^ = 148 (dotted), 153 (solid), and 156 (dashed) GeV, where 
we have fixed the parameters as 

A3 = -1.26 , A4 = -0.0205 , 73 = 11.3 , (35) 

m^± = m^o = m^o^ + 4 GeV , = 135 GeV , rrih = 125 GeV . (36) 

As we see from Fig. \5\ the effective cross section around the decoupling temperature x ~ 20 
has a correct size and decreases drastically at low temperature. The effective cross section is 
normalized to 10~^^ cm^s~^, because it is the size to obtain the observed relic density of DM. 
The scalar couphngs X2,1i and 74 do not enter into the cross sections (ISTl) . and 75,75,77 
and n are irrelevant because is much heavier than rj and x- 

In Fig. [6] we show the area in the 72-73 plane in which the total relic den- 
sity Qrh^ = 0.1157 ± 0.0046 {2a) with = 135 GeV can be obtained for 
m^o = 148 (black circles), 152 (red diamonds), 153 (red circles), 153.3 (blue circles), 
153.5 (green squares), 154 (violet crosses) and 156 (brown triangles) GeV, respectively. The 
right-handed neutrino masses, M^, are all 1 TeV, and the Yukawa couplings are chosen to 
yield J2ik l^ifcP = (10""^)^. The vertical (black dashed) lines are the upper bounds of 72 set 
by the XENON 100 (right) and IT (left) experiments fl2^ . From Fig. E] we can also see that 
there exits a parameter space with = 135 GeV, m^o > 153 GeV, a large 73/47r (between 
0.65 and 1.0) and 72 satisfying the XENON constraint fl29|) . A large 73 is needed to explain 
the 135 GeV 7-ray line in the Fermi spectrum, as we will see below. 



^ The mechanism has been also used in the model of 
LAT. But the light charged scalar faces a problem in explaining the neutrino mass. 



40( 1 explaining the monochromatic 7 at the Fermi 
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FIG. 5: The effective anniliilation cross section (tfie expression in [ ] of ([M]) ) as a function oi x = fi/T for 
m^o = 148 (dotted), 153 (solid), and 156(dashed) GeV with fixed at 135 GeV. 




FIG. 6: The area with nrh^ = 0.1157 ± 0.0046 {2a) and = 135 GeV in the 72-73 plane for 
m^o^ — 148 (black circles), 152 (red diamonds), 153 (red circles), 153.3 (blue circles), 153.5 (green squares), 
154 (violet crosses) and 156 (brown triangles) GeV, respectively. The vertical (black dashed) lines are the 
upper bounds of 72 set by the XENON 100 (right) and IT (left) experiments (|29)) . 



12 




7>7 



■1, 7 



FIG. 7: One-loop diagrams for xx ~^ 11- 
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FIG. 8: S-channel diagrams for h — ?> 77. 



D. Fermi LAT 135 GeV 7-ray line 



Our model contains the coupling between x and r^"*" (£^1 



-73 '7''''7X^) 



Because of 



this coupling there are diagrams in Fig. [TJ which produce monochromatic 7 lines through the 
annihilation of ;^We would like to use these diagrams to explain the monochromatic 7-ray 

There exist also s-channel diagrams 



31 



32 . 1351 . 1361 ] observed at the Fermi LAT 



27-30 



line 

with the SM Higgs propagator as shown in Fig. [81 The hxx couphng is proportional to 72, 
while hrj^r]^ coupling is proportional to A3. The annihilation cross section o'l'q^rj^ — )■ 77) is 
large, because A3 is large (see and ^. Besides, due to the gauge interactions (see 
the right diagram of Fig. H]) the relic density of rjp, is very small so that the annihilation of the 
rjpi^ DM can not contribute to the monochromatic 7-ray. Furthermore, because of the same 
reason, the tree level annihilations of r/^ into a pair of W^s and Z's, which would contribute 
to the continuum 7-ray spectrum, are also suppressed. In contrast to this case, the pure 
gauge interaction (the right diagram of Fig. H]) is absent for the annihilation of x- The entire 
annihilations of x into the SM particles are controlled by the single coupling 72, which has 
to satisfy the constraint fl2^ . Therefore, we may assume that the main contribution to 



In [371] . (yivRVR ~^ 77) is used to explain the monochromatic 7-ray line, where a fine-tuned cancellation 
mechanism to suppress the total annihilation cross section of 77^ [56| is employed 
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FIG. 9: The area for a{xx ll)'" = (l-27Io:43) x 10"^'^ cm^s"^ in the m^±-73 plane. 



77) comes from the one-loop diagrams in Fig. [7] and find 



73 



47r, 



4.23 
2.66 
1.87 
1.39 



X 10 cm^s ^ for m^± = < 



140 
145 
150 
155 



GeV 



(37) 



These values should be compared with (1.27l!o;|3) ^ 10^^'' cm^s~^ 32j for an Einasto DM 
galactic halo profile, which is the size that could explain the monochromatic 7 line observed 
at the Fermi LAT. So, if 73/471 is large of 0(1), the desired value could be obtained. In Fig.[9] 
we show the area in the m^±-73 plane in which o"(xx ~^ 77)^ = (1.27lo;43) x 10"^'^ cm^s~^ 
can be obtained. If rj^ is lighter than then x can be annihilated into a pair of 77^. For 



135 GeV the annihilation cross section becomes ct(xx ~^ V^V = (73/27rm^)(l 



m^±/m^)^/^ — ^K-'- ~ ""^^i/""^^)^^^ X 10" 
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s —1 

cm s , 



which is too large to obtain a sufficiently 



large relic density of x for a large 73 unless rriy. < m^± . This is why we have to assume that 
< m^o^, m^o, m^±. Comparing Fig. with Fig. 0, we see that there is an overlapped area. 



that is, a parameter space in which a{xx ~^ 77)^^ ^ 10 x 10 cm^ s ^ and Qxh"^ 



be obtained. This is shown in Fig. [TO] for m^o^ 
GeV. 



153.3 GeV and m^± 



m„o +4 GeV 

'IR 



0.12 can 
157.3 



The same diagrams as in Fig. [7] also produce 7^. The annihilation cross section is given 



by 



(^{XX lZ)v 



1^ 3 4 \^Q\^ri±Mx)\ 



(38) 
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FIG. 10: The parameter space (overlapped area) in which a{xx ~^ 77)w ~ (1.27^013) x 10 cm^s ^ and 
nrh'^ = 0.1157 ±0.0046 (2cr) can be obtained for m„o = 153.3 GeV and to„± = m„o +4 GeV = 157.3 GeV. 

^ ' 'IR I 'IR 



where 



1 
2 



1 pI—x f 

dx / " 
JO 



Tz{x'' + y^)+Axy{l-Tz/2) 
+ rz{x^ + 2/2) - Axy{l - tz/2) 



(39) 



with Tr, 



'V ~ "^^±/'"^x ~ m|/m^. For m^± = 150 GeV and = 135 GeV, for instance, 

we obtain a(\Y — >■ 7^)f = 7.9(73/47?)^ x 10"^^ cm^s^\ which is about 4 % of cr(xx ~^ ■yy)v. 
According to 57| this is welcome to explain the Fermi LAT monochromatic 7 line. 



We also have to satisfy the constraints on the continuum 7 and follow the analyses of 



57 



fsee also 



58, 



59 



which give similar constraints). They consider two different constraints, 
supersaturation constraint and shape constraint, which can be transferred to the upper 
bound on the ratio of the total annihilation cross section to cr(xx ~^ 77)- Specifically, they 
consider the constraint on the theoretical ratio 



R 



th 



(40) 



where (Tt is the total annihilation cross section. Since the dominant origins for the continuum 
7 are cr{xx ~^ W~^W~ , ZZ,bb,r^T~ , we also consider the individual cross sections 

and calculate 

th ^ixX a) a{xx ^ «) 



R. 



2cr{xX 77) + crixX ^ 7^) 2a(xx ^ 77) ' 



(41) 



in addition to R^ ^ 



(^ixx SM)/2cr(xx 77), where a 
same parameter values given in fl5Bl) with 72/47r 
we obtain 



1.0 X 10 ^ and m„o 

'iR 



W+W-, ZZ, ff. We use the 
153.3 GeV, and 



0.981 X 10' 



0.1197 



0.1197 



(42) 
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and 







1.2 X 10" 


-27 


cm'^s~ 


-1 








r\~i 


8.0 X 10" 


-29 


cm'^s" 


-1 


pth 


~ 0.06 , 


crixX ^ W+W-)v 




3.9 X 10" 


-29 


cm'^s" 


-1 


pth 


~ 0.03 , 


(^ixx -> ZZ)v 


r\~i 


1.7 X 10" 


-29 


cm^s" 


-1 


pth 


~ 0.01 , 


(^ixX hh)v 




2.5 X 10" 


-29 


cm'^s" 


-1 


pth 


~ 0.02 , 


^{XX ff)v 




1.1 X 10" 


-31 


cm'^s" 


-1 


pth 


- 10~^ 



(43) 
(44) 
(45) 
(46) 
(47) 
(48) 



These values of R should satisfy the supersaturation constraint as well as the shape con- 
straint of [57| (see also 



58 



There are also constraints coming from the anti-proton-to- 
proton flux observed by the PAMELA 60|. Anti-protons can be produced by the DM 
annihilations into the gauge bosons, Higgs bosons and quarks. To explain the PAMELA 
data, these productions have to be suppressed. The values of i?*^ given i n (H^- fHS) sa tisfy 
all the constraints including the most stringent one i?*^ ^ 10~^^ cm'^s"^ j57H59l. I6ll-l64|. So, 
the model could explain the monochromatic 7 line observed at the Fermi LAT if m^± 153 
GeV and 73/4?? ~ 0(1), which is at the border of perturbation theory. 



IV. CONCLUSION AND DISCUSSION 

In this paper we have proposed a non-supersymmetric model of two-loop radiative seesaw, 
in which the lepton number is softly broken by a dimension two operator, and the tree-level 
Dirac mass is forbidden by Z2 x This discrete symmetry can be used to stabilize two or 
three dark matter particles. The model contains, in addition to the SM Higgs field, an inert 
SU{2)l doublet scalar rj and two inert singlet scalars (p and and this is a minimal set. We 
have considered the SM Higgs boson decay into two 7's and found that it is enhanced by 
rj~^ circulating in one-loop diagrams for h — 77. 77"^ is also circulating in similar one-loop 
diagrams contributing to XX ~^ 77; "we have found that the model has a potential to 
explain the Fermi LAT 135 GeV 7-ray line. 

The mechanism to explain the Fermi LAT 135 GeV 7-ray line in the present model is 
strongly based the fact that there are at least two particles, particles (which can be DMs, too) 
and one DM particle of similar masses. Let us briefly outline the mechanism. Annihilation 
(or decay) of DM into 7's happens always at the loop-level. Those into the SM particles, 
i.e. W, Z, Higgs boson pairs etc, are usually possible at the tree-level, and they produce 
continuum 7 rays as well as anti-protons. To explain the 135 GeV 7-ray line, we have to 
suppress these tree-level processes somehow, or enhance the loop process, while keeping 
the relic abundance of DM at the observed value. In the present model this is realized in 
the following way. There are two kinds of the tree-level DM annihilations; one into the 
SM particles and the other into a pair of other DM particles (DM conversion). The slightly 
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heavier DM (77^ in our model) has large annihilation cross sections both into the SM particles 
and other DM particles so that its relic density is very small. The 135 GeV 7-ray line comes 
mainly from the annihilation of the slightly lighter DM {x in our model). Its annihilation 
cross section into the SM particles has to be sufficiently small to suppress the continuum 7 
rays and the production of anti-protons. Although the annihilation of the lighter DM into 
heavier DM is kinematically forbidden at zero temperature, this process becomes operative 
at high temperature: The smaller the mass difference of two DM particles is, the more 
effective is the DM conversion. So, at high temperature the annihilation of the lighter DM 
is controlled by the mass deference and can be large, but at low temperature, this conversion 
process practically disappears (see Fig. |5]). (For this mechanism to work, the slightly heavier 
particle does not have to be stable. ) It is however important that the lighter one is SM 
gauge singlet to avoid that the tree-level annihilations are entirely controlled by the SM 
gauge interactions. (In the present model, the x DM is not ad hoc introduced.) This is 
the reason why we can obtain the observed value of the relic density for the lighter DM, 
although the annihilation cross sections of the lighter DM into the SM particles are very 
small ^ 0(10"^'') cm^s~^ in the galaxy. 

The annihilation cross sections into the SM particles given in fH2|) . ( |44|) - (148|) are obtained 
without one-loop corrections. Strictly speaking, we should include the one-loop corrections, 
because the tree-level contributions are so small that the one-loop corrections may be larger 
than the tree-level corrections. Similarly, the relic densities for Fig. and also fH^ have 
been computed by neglecting the co-annihilation of the DM particles with the charged and 
CP odd components of 1], although we have assumed that their mass differences are not 
large (i.e. m^± — = (17 — 25) GeV, m^± — m^o = 4 GeV and m^± = m^o). If the one- 
loop corrections would change 72 effectively by an order of magnitude, then the annihilation 
cross sections ( 1i4|) - fj48|) would change by two orders of magnitude. Even in this case, the 
supersaturation constraint and shape constraint would be satisfied. The co-annihilations 
also would effectively increase 73. To obtain a realistic relic abundance for the x DM in that 
situation, its mass should be slightly increased as one can see from Fig. [HI The one-loop 
analysis including the co-annihilations is beyond the scope of the present paper, and we will 
leave it for our future project. 

We have assumed throughout that the is so heavy that it decays into a 77^ and a x- 
In the case that m^^ < + m^o , it becomes the third DM, whose annihilation may be 
responsible for the second monochromatic 7-ray line in the Fermi data [34]. We leave this 
question to the future program. 
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